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Linear waves in bounded inviscid fluids do not generally form normal modes with regular 
eigenfunctions. Examples are provided by inertial waves in a rotating fluid contained 
in a spherical annulus, and internal gravity waves in a stratified fiuid contained in a 
tank with a non- rectangular cross-section. For wave frequencies in the ranges of interest, 
the inviscid linearized equations are spatially hyperbolic and their characteristic rays 
are typically focused on to wave attractors. When these systems experience periodic 
forcing, for example of tidal origin, the response of the fluid can become localized in the 
neighbourhood of a wave attractor. In this paper I define a prototypical problem of this 
form and construct analytically the long-term response to a periodic body force in the 
asymptotic limit of small viscosity. The vorticity of the fluid is localized in a detached 
shear layer close to the wave attractor in such a way that the total rate of dissipation 
of energy is asymptotically independent of the viscosity. I further demonstrate that the 
same asymptotic dissipation rate is obtained if a non- viscous damping force is substituted 
for the Navier-Stokes viscosity. I discuss the application of these results to the problem 
of tidal forcing in giant planets and stars, where the excitation and dissipation of inertial 
waves may make a dominant, or at least important, contribution to the orbital and spin 
evolution. 



1. Introduction 

Geophysics and astrophysics give rise to problems involving waves in axisymmetric 
rotating fluid bodies. The fluid outer core of the Earth supports inertial waves, for which 
the Coriolis effect provides the restoring force. Such waves have been tentatively identi- 
fied in gravimetric data and related to theoretical work as well as laboratory experiments 
(e.g. Aldridge & Lumb 1987). HeHoseismology uses the observed frequencies of solar os- 
cillations to deduce the internal structure and rotation of the Sun, while asteroseismology 
applies related techniques to distant stars (e.g. Christensen-Dalsgaard 2002). Waves in 
astrophysical accretion discs have been studied because of their important role in the 
interaction of a planet or other satellite with the disc in which it forms (e.g. Goldreich 
& Tremaine 1980), and also in an attempt to explain observed quasi-periodic oscillations 
from accreting white dwarfs, neutron stars and black holes (e.g. Kato 2001). Some of 
these problems involve free modes of oscillation that may grow through mechanisms of 
overstability or may be excited by turbulent noise. Tidal interactions between orbiting 
and spinning bodies, however, involve an almost strictly periodic forcing of waves. 

After the linearized wave equations governing the dynamics of the fluid are reduced 
by Fourier analysis in time and azimuth, there results a two-dimensional problem for the 
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spatial structure of the wave in the meridional plane (see Appendix El- In the absence 
of dissipative effects such as viscosity, this problem may have either elliptic or hyperbolic 
character depending on the frequency of the wave. A mixed-type problem may also 
occur if the character of the equations changes from one part of the body to another. 
Typically, hyperbolic character occurs when the wave frequency (Doppler-shifted, in the 
case of a non-axisymmetric wave, into the frame rotating with the local fluid) is small, 
in the range of inertial and gravity waves rather than acoustic waves. When a hyperbolic 
system of equations is posed in a finite region with physical boundary conditions on a 
closed surface, rather than the mathematically canonical Cauchy boundary conditions on 
an open surface, the problem can be regarded as ill-posed and generally does not possess 
a regular solution. This property presents a serious difhculty when one seeks the response 
of a body to low-frequency tidal forcing. f 

It is possible to sidestep this issue in certain problems of special symmetry where a 
separation of variables can be applied. For example, in a spherically symmetric non- 
rotating star the linearized equations can be projected on to spherical harmonics and the 
problem reduced to a system of ordinary differential equations in the radial direction. 
Most treatments of tides in stars and giant planets have used this approach (e.g. Zahn 
1970), or have included the Coriolis force within the so-called 'traditional approximation' 
(in which only the radial component of the angular velocity is considered), which also 
permits a separation of variables (e.g. loannou & Lindzen 1993; Savonije & Witte 2002). 
The simplest problems involving inertial waves, featuring an incompressible fluid in a full 
spherical, spheroidal or cylindrical container, can also be solved by separation of variables 
and appear to possess complete sets of normal modes with regular eigenfunctions, even 
in the absence of viscosity (e.g. Greenspan 1968; Zhang, Liao & Earnshaw 2004). 

Unfortunately the approximation of a non-rotating star is rarely valid in problems 
of tidal forcing. The tidal frequencies are linear combinations of the orbital and spin 
frequencies with small integer coefficients, apart from small corrections due to any pre- 
cessional effects. In most problems of interest, the frequencies of the most important tidal 
components are smaller in magnitude than twice the spin frequency of the body, and the 
Coriolis force cannot be neglected (Ogilvie & Lin 2004). The traditional approximation 
is valid only in stably stratified regions where the tidal frequency is much less than the 
buoyancy frequency, and does not apply in convective regions of stars or giant planets 
where the buoyancy frequency is essentially zero. In such regions the wave equations are 
hyperbolic in character, describing pure inertial waves, and cannot be solved by sepa- 
ration of variables. Numerical solutions of the two-dimensional problem were obtained 
by Savonije, Papaloizou & Alberts (1995) and Savonije & Papaloizou (1997) in the case 
of an early-type star with a small convective core and, more recently, by Ogilvie & Lin 
(2004) in the case of a giant planet with an extended convective region. 

Numerical studies have revealed the intricate structure of inertial waves in an incom- 
pressible rotating fluid contained in a spherical annulus (HoUcrbach & Kerswell 1995; 
Rieutord & Valdettaro 1997). In order to find normal modes it is necessary to include a 
viscosity so that the problem becomes of elliptic character and is mathematically regu- 
larized. As the viscosity tends to zero the eigenfunctions become increasing localized in 
the neighbourhood of singular linear structures known as wave attractors, as described 
in greater detail by Rieutord, Georgeot & Valdettaro (2001) and Rieutord, Valdettaro & 
Georgeot (2002). A wave attractor can be understood as a limit cycle towards which the 

f In accretion discs, where the systematic shear causes a strong Doppler shift of the frequency 
of non-axisymmetric waves, the equations typically have hyperbolic character in the vicinity of 
the corotation resonance. 
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characteristic rays of the inviscid wave equation are focused as they reflect repeatedly 
from the boundaries of the container. Unhke the case of a full sphere, it appears that 
inertial waves in a spherical annulus do not possess regular cigcnfunctions in the absence 
of viscosity, apart from the exceptional 'r modes' or 'toroidal modes' which involve no 
radial motion. The presence of an inner core introduces a complexity into the reflection 
patterns of the characteristic rays and causes them to be focused on to wave attractors. 
The existence and importance of these closed ray circuits in spherical shells have been 
known for some time (Stern 1963; Bretherton 1964) and many features of this problem 
were understood before the advent of high-resolution numerical calculations (Stewartson 
& Rickard 1969; Stewartson 1972). 

Wave attractors have also been studied for internal gravity waves. In a notable ex- 
periment, Maas ct al. (1997) applied a periodic forcing to a narrow tank containing 
a stratified salt solution. For wave frequencies smaller in magnitude than the uniform 
buoyancy frequency of the fluid, the characteristic rays of the inviscid wave equation 
are straight lines with a definite angle of inclination depending on the frequency. In an 
upright rectangular tank the rays would propagate around the container ergodically or, 
for special frequencies related to the rational numbers, would close after a finite number 
of reflections. By making one of the side-walls sloping, Maas et al. (1997) introduced a 
wave attractor, in this case an inclined rectangle, into the problem (Maas & Lam 1995). 
Their experiment indicates that the forced response is localized in the neighbourhood of 
the attractor, and this concentration ultimately leads to secondary effects such as mixing. 
Inertial waves in a rotating tank of similar shape have also been studied by Maas (2001) 
and Manders fc Maas (2003). 

There is a close analogy between these problems. The relation between the tank with 
the sloping side-wall and the upright rectangular tank is similar to the relation between 
the spherical annulus and the full sphere. In each case rays in the geometrically simpler 
container (rectangular tank or full sphere) propagate ergodically or, for special frequen- 
cies, form periodic orbits. The inviscid wave equation admits regular eigenfunctions for 
a countable set of frequencies, which may or may not be those for which the rays form 
periodic orbits. In the geometrically more complicated, but more generic, containers the 
rays are focused on to one or more wave attractor for almost all wave frequencies, and 
the inviscid wave equation does not admit regular eigenfunctions. 

In problems of tidal forcing it is of great interest to know how the total dissipation 
rate varies with the forcing frequency, as this determines the rate of secular evolution 
of tidally interacting systems. When the viscosity is small, systems that possess regular 
inviscid normal modes can be expected to exhibit a strong resonant amplification of 
the dissipation rate in the vicinity of the eigenfrequencies, but very small dissipation 
elsewhere. Although the eigenfrequencies may be everywhere dense in some interval, 
the smooth forcing will have a significant overlap only with a few of the lowest-order 
modes. In contrast, wave attractors have a structural stability and exist in intervals of 
frequency that depend on the geometry rather than the viscosity. Systems possessing 
wave attractors can be expected to (exhibit a richer response with significant dissipation 
occurring over extended ranges of frequency. 

The purpose of this paper is to analyse the linear response of such systems to periodic 
forcing in the low-viscosity limit relevant to the geophysical and astrophysical applica- 
tions. I work with simplified model problems but also explain how they are related to the 
original systems. The main result is a validation of the conjecture of Ogilvie & Lin (2004) 
that, in problems where a wave attractor occurs, the total dissipation rate tends to a 
non-zero value that is independent of both the magnitude and the form of the small-scale 
damping process of the waves. This outcome diff'ers significantly from problems that pos- 
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sess regular inviscid normal modes and is of primary importance for tidally interacting 
systems. The analysis presented here provides a method of calculating the asymptotic 
dissipation rate and also describes the spatial form of the forced disturbance. 

The remainder of this paper is organized as follows. In ^211 briefly describe the problem 
of forced internal gravity waves as a definite example of the system under consideration. 
I define a prototypical problem in ^ and construct the asymptotic solution in the limit 
of small viscosity. In f^ll adapt the analysis to a related problem in which the damping 
mechanism is not of a viscous nature. I present the results of direct numerical calculations 
in inland compare these with the asymptotic theory. A concluding discussion is found in 

m 



2. Forced internal gravity waves 

In this section I outline the preliminary analysis for forced internal gravity waves in a 
narrow tank. This is essentially identical to the problem studied by Maas et al. (1997) 
except that a vortical body force is assumed rather than a parametric forcing. 

Consider a fluid initially at rest in a gravitational field —g and with a uniform 
temperature gradient l3ez, where (x,y,z) are Cartesian coordinates. In the Boussinesq 
approximation the linearized equations for the velocity u and the temperature pertur- 
bation are (Chandrasekhar 1961) 

— = -Vvj + gaO -f i/V^w + a, (2.1) 

f)9 

— = ^(3u, + nVH, (2.2) 
at 

V-M = 0, (2.3) 

where w is the pressure perturbation divided by the reference density, a is the coefficient 
of expansion, v is the kinematic viscosity and k, is the thermal diffusivity. The fluid is 
excited by an external body force a per unit mass. 

Consider a two-dimensional problem, approximately representative of the situation in 
a narrow tank, in which Uy — ay = and all quantities are independent of y. The velocity 
is described by a streamfunction ip such that 

dip dip 
= u.= g^- (2.4) 

I focus on the limit of large Prandtl number in which k can be neglected (except in thermal 
boundary layers, which are of negligible thickness compared to the viscous boundary 
layers). In reality, the Prandtl number of water exceeds 10 for temperatures below about 
10°C, so this limit may be of some relevance for laboratory experiments or for terrestrial 
tides confined in oceanic basins. 

When considering the long-term response to a periodic force with (real) angular fre- 
quency w, all perturbation quantities may be assumed to have the form 

M Re [M(a;, z) e"'"'] , (2.5) 

etc. Eliminating tSj and 0, one obtains the linearized vorticity equation in the form 

i^V2^ = — ^-^V^V^ + Z, (2.6) 
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where N'^ = ga(3 is the square of the buoyancy frequency and 

/ = f^-f^ (2.7) 
az ox 

is the y-component (and only non-vanishing component) of the vorticity forcingf V x a. 
With rigid waUs the boundary conditions are 

= 0, n • = 0, (2.8) 

where n is the outward normal vector on the boundary. The time-averaged dissipation 
rate (per unit length in the y-direction) is 

D = ]^jjv\V'^^\^dxAz, (2.9) 

where the integral extends over the area of the container in the a;z-plane. Using equation 
H2.5|) and the boundary conditions one can relate the dissipation to an overlap integral 
between the vorticity forcing and the streamfunction in the form 

D =^Re JJ i;*fdxdz. (2.10) 

The spectrum of internal gravity waves lies in the interval —N < to < N. In this range 
of frequencies the inviscid version of equation H2.6|l . with 1^ = 0, involves a hyperbolic 
operator 

The characteristics of the inviscid equation are straight lines with slopes ±aj(-/V^— w^) 
that depend on the wave frequency. The occurrence of wave attractors in containers of 
various shapes has been illustrated by Maas et al. (1997). 



3. A prototypical problem and its asymptotic solution 

3.1. Definition of the problem 

In this section I slightly abstract the above problem and define a prototypical model 
problem that is very closely related but marginally simpler to solve. Consider the equation 

for an unknown function ^p{x,y) in a domain V C R^, with boundary conditions 

iP = 0, n • V-0 on dV. (3.2) 

Here e ^ 1 is a small dimensionless parameter, /(x, y) is a given complex-valued function 
and n is the outward normal vector on dV. 

As suggested by ^ one can think of ip as the streamfunction of a two-dimensional flow 
of an incompressible fluid subject to no-slip boundary conditions. The viscosity of the 
fluid is proportional to and the curl of the body force is proportional to /. The limit 
e 1 corresponds to the Reynolds number being large. The problem is linearized and a 
harmonic time-dependence of all quantities has been assumed. The fluid has a restoring 

f It may be noted that tidal forcing always derives from a potential and therefore has no curl. 
Nevertheless, the dynamical tide in a body of non-uniform density, or one with a free surface, 
does experience a vortical effective forcing as explained in Appendix IbI 
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Figure 1. Example of a quadrilateral domain containing a unique wave attractor. 

force such as buoyancy or the Coriolis force. In fact the inviscid part of the problem of 
forced internal gravity waves considered above can be rendered exactly in this form by 
a linear transformation of the coordinates, depending on the wave frequency, that maps 
the sloping characteristics on to horizontal and vertical lines. In the process the viscous 
operator is slightly modified, although this has essentially no practical consequences. 

Consider the inviscid version of equation (|3.1I) . obtained by setting e = and aban- 
doning the no-slip boundary condition n • V?/' — 0. The inviscid equation is hyperbolic 
and its characteristics are horizontal and vertical lines. In general, the inviscid problem 
is mathematically ill-posed and does not possess a regular solution. Instead, one must 
consider the limit e — > of the viscous problem. 

The time-averaged energy dissipation rate may be defined as 



Using equation (|3.1() and the boundary conditions (|3.2|) one can relate the dissipation to 
an overlap integral between the vorticity forcing and the streamfunction in the form 



The problem at hand is to determine how D depends on e in the limit e — > (and also 
on / and V if these are varied) . 



A ray segment consists of a horizontal or vertical line segment connecting two points 
of dT) {vertices) and lying wholly in V. A ray circuit consists of a connected set of 
consecutive ray segments and represents the propagation of a wave characteristic around 
the container. Typically a ray circuit will be of infinite length. Exceptionally, a circuit 
that closes on itself after a finite number of reflections is a periodic orbit. 

In principle, the domain V may contain any number of periodic orbits. I consider 
the case in which 7) contains exactly one periodic orbit and the orbit is simple (i.e. a 
rectangle). Without loss of generality one may place the origin of coordinates at the 
bottom left-hand corner of the rectangle so that the corners OABC are at (0, 0), (X, 0), 
{X,Y) and (0, F) with X,Y >Q.l label the segments 1, 2, 3, 4 in the order OA, AB, 




(3.3) 




(3.4) 



3.2. Ray circuits and the wave attractor 
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BC, CO. Let the slopes of the boundary at the four corners be = ^i, sb = , 
sc = ts and so = —t^^. Then tj > is the tangent of the angle between the 'incident' 
ray (in the positive sense OABC) and the boundary at the end of segment j. 

For rays close to the periodic orbit, the boundary may be represented locally by linear 
approximations 

y ~ ti{x — X) near A, 

y w ~t^^(x - X) + Y near B, 

y w t^x + Y near C, 

y ~ —t^^x near O. (3.5) 

The local approximate mapping of vertices defined by the ray propagation (in the positive 
sense OABC) is then 

{-t4y,y) ^ {ti^y + x,y) 

^ {t^^y + X,-t^h^_^y + Y) 

^ i-t^\%'y,-t^\'y + Y) 

^ i-tl%%'y,t^%%'t^'y)- (3.6) 

On each loop, therefore, 

{-tiy,y) ^ a~^{-t4y,y), a = tit2tzti. (3.7) 

It may be expected that a ^ 1 in a typical asymmetric container. The periodic orbit is 
then attracting, and focuses rays propagating either in a positive sense (if a > 1) or in 
a negative sense (if < a < 1). I assume without loss of generality that a > 1 (the case 
< a < 1 can be obtained by a reflection of the problem). The focusing direction is then 
OABC and a is the focusing power of the attractor. 

A particularly useful illustrative example is provided by a quadrilateral container (fig- 
ure P). The quadrilateral is defined uniquely by specifying the rectangle OABC and the 
four tangents tj . It can be shown that the wave attractor is unique. All ray circuits other 
than the attractor itself converge towards the attractor at either end. I assume that this 
property holds in the domain P. 

The attractor divides the boundary into two zones, one consisting of OEA together 
with BCC and the other consisting of AFB together with CHO. Ray circuits other than 
the attractor itself have all their vertices in one zone or the other. Therefore the attractor 
has two 'sides' to it, and any ray circuit converges towards the attractor from the same 
side at either end. 

For rays near the attractor one can measure the distance from the attractor using the 
y-coordinates of the appropriate vertices (the ones near O). The entire bundle of rays 
associated with one side of the attractor (e.g. the 'positive' side y > 0) can be labelled 
using a variable 9^ G [0, 1) defined as follows. One starts a ray at a vertex very close to 
O and with y-coordinate y — a^"^^ for some arbitrarily large positive integer n » 1. 
The ray is followed in the negative sense as it spirals away from the attractor. Eventually 
it turns around and converges towards the attractor again, having vertices very close to 
O with y-coordinates y ~ a~"~® n 3> 1. This mapping between one end of a ray 

circuit and the other defines a function 0+(6'+), with 8"^ S [0, 1). 

The variable 0^ is periodic in nature (6'"'" = 1 is identified with 6'+ — 0) since the 
ray bundle wraps around itself. The function 9+(0+) can be considered topologically as 
a map of the circle on to itself. Owing to the reversibility of the ray circuit it has the 
self-inverting property 9+(9+(0"'")) = which will be used later. 
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For the negative side y < of the attractor there is a similar variable 0~ and a mapping 
Q~ {6~). In problems with multiple or non-simple attractors a more complicated mapping 
of ray bundles is likely to exist. 

3.3. Construction of an asymptotic solution 

In the limit e — > the solution divides into three regions: (i) an inner region localized near 
the wave attractor and consisting of a detached shear layer of width 0(e); (ii) an outer 
region consisting of most of the rest of the container; (iii) standard viscous boundary 
layers of width 0(e^/^) close to the walls. There are also boundary layers in the corner 
regions close to the vertices of the attractor. 

I first simplify the problem by assuming that a particular solution for the inviscid 
problem can be found, i.e. a function iJj{x, y) satisfying the equation 

without regard to boundary conditions. (This can be obtained by an indefinite integration 
of — i/ with respect to x and y.) The desired solution is then %p = ip -\-jp where ip satisfies 
the equation 



dxdy 

and the boundary conditions 



+ e'^VV = -e-'V^V (3.9) 



■tp^-i^, n-V-p ^ -n-Viij. (3.10) 

Note that / and ip are supposed to be smooth and do not have any fine structure associ- 
ated with boundary layers or shear layers. Therefore the right-hand side of equation 1)3. 9|l 
is small everywhere and will be unimportant in constructing the leading-order asymp- 
totic solution. The inhomogeneity has effectively been transferred from the differential 
equation to the boundary conditions. 

3.4. The outer solution 
The asymptotic outer solution is simply of the form 

^^^(°"*)(x,2/). (3.11) 



It satisfies the inviscid problem 



with the boundary condition 



a2 T(out) 



The second boundary condition is taken care of by the intervention of a standard viscous 
boundary layer (region (iii) mentioned above). 
The general solution of equation H3.12|l is 

=5(2;)-%), (3.14) 

where g and h are functions to be determined. Consider two consecutive vertices, P and 
Q, say, of a ray circuit. If they are connected by a horizontal ray segment, they share the 
same value of y and therefore of h{y). It follows from equations H3.13|l and H3.14|l that 
gixg) — g(xp) — tpp — tpQ. Similarly, for vertices connected by a vertical ray segment, 
^(uq) ~ h{yp) = ipQ — "pp- Since ■)/' is known, a knowledge of g (or h) at any vertex would 
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be sufficient to propagate the solution along the ray circuit. However the given boundary 
conditions do not provide such information. 

Consider the ray circuit as described in i )3.2l At each end of the circuit the ray indefi- 
nitely repeats a loop around the attractor OABC. Around each loop the value of g must 
change by an amount 

S = ipo - 1pA + ipB - ipc, (3.15) 

and the value of h also changes by S. Since ^ is smooth, one need not worry about the 
small variation of tp close to the vertices of the attractor. In fact 

6 = [ [ ^-^dxdy = -i [ fdA, (3.16) 



dxdy 

where the integration is over the area enclosed by the attractor. Therefore 6 is indepen- 
dent of the choice of particular solution ^. It is also notable that if, as in §21 / is the 
perpendicular component of the curl of the external force per unit mass, then Stokes's 
theorem provides a simple relation between S and the line integral f a • dr of the force 
around the attractor. 

It follows that h{y) exhibits an increasingly fine structure and, in general, a logarithmic 
divergence as y — > 0. The same is true of g{x) near x = and x = X, and of h{y) near 
y = Y. This singular behaviour of the outer solution is the reason that an inner region, 
in which viscosity regularizes the solution, is required. For small y, one has 

h{a-^y) - h{y) K S. (3.17) 

Consider the positive side y > first and define the variable y^ = — Iny/lna, which 
increases by 1 on each loop. Then h{y) ~ h'^{y'^) as ?/ \ 0, with 

h+{y+ + l)-h+{y+)^5. (3.18) 

The general solution of this functional equation is 

h+{y+)=y+5 + H+{e+), (3.19) 



I.e. 



where 



h(y)^ \-p^]5 + H+{e+) asy\0, (3.20) 



In 



0+ = y+ n^od I = (^-^^ mod 1 (3.21) 

is the ray-labelling variable introduced previously and is an arbitrary function with 
period 1. It may be expanded in a Fourier series 

oo 

H+{0+) ^ J2 H+e^"'''''*. (3.22) 

Tl — — CX3 



hiy)-^[-p^]d + H-ie-) as2//0 (3.23) 



For the negative side of the attractor, a similar result holds, with 

In \y\ 

In a 
and 

oo 

H-{9-)= iJ-e^"^'^". (3.24) 
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Now consider the total change in h(y) around the ray circuit, from a vertex with 
y = a~"i^^ , rii ^ 1, to a vertex with y ~ a~"^~® \ n2 ^ 1- This change can 
be computed numerically by summing the values of tJj, with appropriate signs, at the 
vertices visited by the ray circuit. The total change can be written as 

[ns + e+{0+) ^ni-e+]S + A+{e+), (3.25) 

where A"*" does not depend on ni or n2 because to increment either by 1 adds a loop 
around the attractor at one end and changes h(y) by at that end. Owing to the 
reversibility of the ray circuit, A+(e+(6'+)) = -A+(6'+). 

Connecting the asymptotic forms, equation (|3.2Q|I . of the solution at the two ends of 
the circuit, one finds 

H+{e+{9+))-H+{9+) = A+{9+). (3.26) 
The general solution of this functional equation is 

H+ie+) = -^A+{0+) + J+{0+), (3.27) 

where J"*" is an arbitrary function satisfying the symmetry condition J+(0+(0+)) = 
J~^{d~^). An analogous equation holds for the negative side of the attractor. There is 
insufhcient information to determine the solution uniquely until the equations for 
and can be linked through the dynamics in the inner region. 



3.5. The inner solution 

The inner solution is localized near the attractor and consists of four segments. I begin 
by defining coordinates parallel (^) and perpendicular (77) to the ray in each segment: 



6 




m 






(sej 


jment OA), 


6 


= y, 




= e-\x- 




(sej 


jment AB), 


6 


= (X - x), 


m 




-y) 


(sej 


jment BC), 


^4 


= Y~y, 


V4 


= e-\-x] 


1 


(sej 


jment CO). 



(3.28) 

In segment j the parallel coordinate runs from to Sj, where Si = S3 = X and 
S2 = S4 = Y . The perpendicular coordinate -qj is stretched to resolve the inner region 
of width 0(e). In segment j the asymptotic inner solution is of the form 



(3.29) 



and satisfies the equation 



which follows from equation 1)3. 9|l under the assumed scalings. 

Part of the difficulty in solving equation (|3.3U|) consists of connecting the four segments. 
Boundary conditions apply in the corner regions where the segments meet and the wave is 
reflected. Where segment j overlaps with segment j + 1, the asymptotic solution consists 
simply of the linear superposition of V'j™'' a-nd V'j+ii because equation H3.30|l is linear 
and homogeneous. In the corner region, and ^j+i are essentially constant while rjj and 
77j+i vary by 0(1). Therefore the corner solution is of the form 

^~V;r(S„r;,)+^(V"l(0.'?.+i)- (3-31) 
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The location of the boundary is given by 

%+i -i7'% +0(e). (3.32) 

From the boundary condition ^ = —ijj one obtains the relation 

ViSl(0,t7S,) = -V^f ^(S„r;,) (3.33) 

which connects the solutions in consecutive segments. (To satisfy the second boundary 
condition n • ^ij) = a thinner boundary layer must intervene in the corner.) At the 
vertex O the corresponding relation is 

^^^(0,t4 ir;4) = -^f ^(S4,r;4) - V'4. (3.34) 

I now define a connected inner solution by concatenating the solutions in the various 
segments. To connect the perpendicular coordinates smoothly, let 

{rji in segment OA, 

tii]2 in segment AB, , , 

^i^2'73 in segment BC, 

tit2t3ri4 in segment CO, 

i.e. rj = fjTjj in segment j, with /i = 1 and 

fj^Yltk, J>1. (3.36) 

k=l 

To leave equation H3.30|l in the same form, the parallel coordinate must be rescaled 
by a factor For a continuous concatenation of the parallel coordinates, let 

^1 in segment OA, 

t _ I /l^ + Si in segment AB, 

^ VI6 + S1+/IS2 ^ in segment i?C, ^^'^'^ 
+ =1 + /|22 + /|=3 in segment CO, 
i.e. ^ = ^1 in segment 1 and 

^ = /|?.+E/'^^ (3.38) 

fc=i 

in segment j > I. This coordinate runs from to S around the attractor, where 

4 

- = E = (1 + + + ^2^3)^- (3-39) 

The solution itself, away from the corners, may then be written in the form 

^~V;(-)(^,,7) (3.40) 

with 

V^^'-He,^) - (-l)^"V?"^(e„^,) +E(-l)'^fe (3.41) 

k=l 

in segment j. The boundary condition (|3.33|) then translates simply into the condition 
that '7) be continuous at each vertex. Meanwhile the equation (|3.3Uf) becomes 

simply 



1- 



+ = 0. (3.42) 



d^dr] drj 
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This concatenation procedure effectively 'irons out' the corners in the attractor, allowing 
the solution to proceed continuously, ft is important to note that the coordinates 77) 
cover the corner regions twice. As noted above, in the corner regions the asymptotic 
solution is the sum of the solutions in the segments that meet there. 

In fact the boundary condition at vertex O is different because the perpendicular 
coordinate rj undergoes a net focusing by a factor of a after one loop around the attractor. 
The boundary condition there is 

V^^'") (S, ar]) = (0, 77) + 6. (3.43) 

To satisfy this condition requires a kind of self-similar viscous spreading of the shear 
layer to compensate for the geometrical focusing. Such a self-similar expansion occurs 
in the solutions of Moore & Saffman (1969) for detached shear layers in rotating fluids, 
which were used by Rieutord et al. (2001) in their analysis of inertial waves in a spherical 
annulus. 

Define a similarity variable 



T 



= (3-44) 

where fi = ^ + c with c to be determined. For r to map continuously at the vertex O we 
require 

{E + c)-'^^^ari = c'^/^ri (3.45) 

and therefore 

(3.46) 



a3 - 1 
One then writes 

^('"Hf,r7) = |<5 + *(M,r) (3.47) 

to accommodate the net increase around the loop. Here /i runs from c to c + S around 
the attractor. Equation H3.42|l translates into 

and the boundary condition (|3.43|) requires simple continuity: 

^'(S + c,r) ^'(c,t). (3.49) 
Making the further transformation 

fi = ce^, (3.50) 
where A runs from to 3 In a, one obtains 

. i d-^ IT 



dXdr 3 ar 3 dr^ dr^ 



0. (3.51) 



The variable A is now periodic and ignorable, so the solutions are found by separation of 
variables to be of the form 

*n = e"="\„(T), (3.52) 
where fc„ — 2n7r/(31na), n G Z and Xn satisfies the ordinary differential equation 



fcn + 3 j Xn - jXn + x'n = 0- (3.53) 

Solutions that do not diverge superexponentially as \t\ — > 00 can be obtained by the 



Wave attractors and tidal disturbances 13 

Laplace transform method as in Moore & Saffman (1969), leading to the integral repre- 
sentation 

poo 

x'^{t)^i e-'P^e-p'p-^'''"dp. (3.54) 



(These functions are related to Moore-Saffman functions of complex order and can be 
expressed in terms of generalized hypergeometric functions.) The asymptotic behaviour 
is 

Xnir) (-3ifc„)! e-(3/2)-'="r-i+3''=" as r ^ +oo, (3.55) 
x'nir) - -(-3ifc„)! e(3/2)-fe" |^|-i+3ifc. as T ^ -oo. (3.56) 

Therefore Xo(''") diverges logarithmically as r — > ±oo, while Xn{T) for n ^ is bounded 
but oscillates indefinitely. One may define XniT) uniquely such that Xn(0) = 0, say. Then 

Xn{r) - lz^eT(3/2).fc„|,|3i.„ _^,± as r ^ ±oo (3.57) 

for some constants c^. In the case n = 0, however, 

Xo(''') In |t| + Cg as T ^ ±oo. (3.58) 
The general solution of the inner problem is therefore given by equation (|3.47|l with 

oc 

^i^l,T)^ A^''"KXn(T)+6„], (3.59) 

n— — oo 

where a„ and 6„ are undetermined coefficients. Noting that |Tp/i = jyyp one finds the 
outer limit of the inner solution to be 



r'''> (^, 7?) -^S + ao In \tj\ - — ln(e + c) + oqCo + bo 



E 



(3.60) 



as ?7 ^ ±oo. 



3.6. Asymptotic matching 

The aim is now to determine the coefficients a„ by matching the outer limit of the 
inner solution to the inner limit of the outer solution. (The coefficients 6„ will not be 
required for a calculation of the asymptotic dissipation rate.) I consider the matching on 
segment 1; it is straightforward to show that the solution then matches in the same way 
on the other segments. 

The inner limit of the outer solution near segment 1 is, for y > or y < respectively, 

^(°-^){x,y) ^ g{Q) - (^-^^y - H^{9^), (3.61) 

where again 

i^^modl. (3.62) 

Ina / 

To compare this solution with the outer limit of the inner solution, equation (|3.60|l . note 
that 77 = e^^y and therefore 

|^|3ifcTi ^2n7vi\n \rj\/ In a g— 2n7ri^^ g — 2n7ri In e/ In a (3 63) 
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I consider first the positive side y > of the attractor and identify 

ao = , (3.64) 

ma 



(-3ifc„)! 

Matching on the negative side yields the same expression for ao, and 

3ikn 



-3ifc„)! 



g-(3/2)^fc„g2n^ilne/lna^-^^^ ^ ^ ^g gg) 



It follows that the functions _ff + (9^ ) and H [6 ) are related to each other in the Fourier 
domain byf 

i7+ = e3-''^"H-, n^O. (3.67) 

Physically, this relation is achieved by a viscous diffusion across the shear layer, which 
allows information to be transmitted across the characteristics of the inviscid equation. 
Now using the identity 



smh(7ra;) 



for real x, one obtains 



|a_„p = ^(|i/+p-|i7-p), n^Q. (3.69) 

3.7. The dissipation rate 

The energy dissipation rate D is of order unity in the limit e and is dominated by the 
inner region. To see this, consider the following order-of-magnitude estimates, bearing in 
mind that the dissipation rate is the area integral of the viscosity times the square of 
the vorticity. In the outer region the streamfunction, velocity and vorticity are all 0(1), 
the area of the region is 0(1) and the dissipation rate is therefore 0{e^); it is small just 
because the viscosity is small. In the boundary layers the velocity is 0(1), the vorticity 
is 0(e~^/^), the area is 0{t^^^) and the dissipation rate is therefore 0(e^/^). Finally, 
in the inner region the streamfunction is 0(1), the velocity is 0(e~^), the vorticity is 
0(e~^), the area is 0(e) and the dissipation rate is therefore 0(1). I now calculate an 
exact asymptotic expression for the dissipation rate. 
The element of area in segment j of the inner region is 

AA^dxAy^e d^^ dr/^ = eJ/M^ dry. (3.70) 

Under the transformation from (^,77) to (/i,T) this becomes 

d^ = e/r4^i/3d^dr. (3.71) 

The vorticity is dominated by perpendicular derivatives of the streamfunction: 

f The mean values of , H^^, are not determined by the matching procedure. The relation 
KUfl impi ies that, if H are analytic functions, they are identical except that their arguments 
are shifted by i7r/ in a in the complex plane and they also differ in value by an additive constant. 
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The dissipation rate at leading order therefore simphfies to 

dAdr. 





/•oo p3\na 




2. 


'-00 Jo 





The different 'modes' proportional to e''"'"'*' add in quadrature with the result 

3 



D 



Ina ^ |a„prfn, 



n— — oo 



where 



dn = / \Xnir)f dr 
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(3.73) 



(3.74) 



(3.75) 



is just a property of the basis functions. The integrals dn are convergent; the dissipation 
is well localized even if the streamfunction is not. To evaluate dn I use the integral 
representation (|3.54l) to obtain 



dr, 



oc fOC />oo / \ 3i/Cn 



-ooJO Jo 



pq dp dq dr. 



-2p „2 



The integration with respect to r is carried out first, using the identity 

e'(9-P)^dr = 2TTS{q-p), 

and leading to the simple result 

dn = 2tt 

The asymptotic dissipation rate is therefore 

oo 

D - \na ^ \anf. 

n— — oo 

Note that |a„p can be related to the outer solution by equation (|3.64f) or H3.69|l . 



p dp 



(3.76) 



(3.77) 



(3.78) 



(3.79) 



4. A problem with a non-viscous damping force 

4.1. Motivation 

In this section I investigate a related problem in which the weak viscous force, which 
depends on second derivatives of the velocity, is replaced by a weak damping force of non- 
viscous form, depending on the velocity itself. I show below that the same asymptotic 
dissipation rate is obtained in this case, implying that the dissipation rate is robust, 
being independent of both the magnitude and the form of the damping force in the limit 
that the damping is weak. 

One reason for investigating this robustness is that the true mechanism by which tidal 
disturbances in astrophysical bodies, such as inertial waves in the convective regions of 
giant planets, are damped is uncertain. The true viscosity is exceedingly small, and the 
waves may decay more readily by other means. Their interaction with convective motions 
might be modelled in terms of an eddy viscosity, but the validity of such an approach is 
unclear and estimates of the relevant effective viscosity vary widely. Alternatively, inertial 
waves may undergo Ohmic damping in the presence of a magnetic field since they couple 
to Alfven waves on small scales. They may also experience nonlinear parametric decay 
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into waves of shorter wavelength. The last two mechanisms are probably modelled more 
accurately, although still only crudely, by using a 'frictional' damping force proportional 
to the velocity, rather than a viscous force. 

Another way to view the following analysis is that it results from the Landau pre- 
scription in which the frequency of the forced wave is given a small positive imaginary 
part, thereby rendering the inviscid problem soluble, and a limit is taken in which the 
imaginary part of the frequency tends to zero. Such an approach might be justified by 
a consideration of the late-time asymptotic solution of the inviscid initial- value problem 
constructed by a Laplace transform method. 

4.2. Sketch of the analysis 
Instead of equation H3.1|l I now consider 

with the boundary condition 

^ = on dV. (4.2) 
The dissipation rate for this problem is 



D^^ l^elV^l^dA^^ReJ^r/dA. (4.3) 



The analysis proceeds in a very similar way to the previous calculation and it is nec- 
essary only to point out the essential differences. One simplification is that no viscous 
boundary layers occur, although in fact it was not necessary to consider the boundary 
layers of the previous problem in any detail. The essential differences appear in the anal- 
ysis of the inner region, which is again of width 0(e). The reduced equation (|3.3()(l in 
segment j becomes 

The factors of and t^ in equations (|3.37t{S35|l should be replaced by single powers 
and the concatenated inner equation (|3.42l) becomes 



52j;(i„) 52j,(i„) 

A similarity solution satisfying the boundary conditions can be obtained in an analogous 
way, now using variables r = r//(^ + c) and A = ln(^/c -I- 1) with c = !B/(a — 1). The 
parallel coordinate A runs from to In a, not 3 In a as previously. The inner equation in 
similarity variables now reads 

and has solutions 

*n = c3'^"^Xn(r) (4.7) 
where I define A;„ — 2n7r/(3 Ina) as previously, and Xn now satisfies 



(4.8) 
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The Laplace transform method leads to 

/>OC 

xUr) = i / e-'P'e-Pp-''''- dp = (-3ifc„)! e-^^/^).^,. _ -y 



l+3ifc„ 
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(4.9) 



where the branch cut is confined to the upper half-plane. The asymptotic behaviour as 
T ^ ±00 is precisely as in equations H3.55|) and H3.56|l . and the asymptotic matching 
therefore yields precisely the same values for the coefficients a„. 
Finally, the asymptotic dissipation rate is computed as 

2 



D 



1 



/"OO pin a 


a* 


/-oo Jo 


dr 



n— — 00 



where now 



dn 



(4.10) 



(4.11) 



The end result is that the asymptotic dissipation rate is exactly the same as for the 
viscous problem, 



D 



(4.12) 



n— — 00 



5. Numerical solutions 

Numerical solutions of equations (|2.6|) . H3.1|l and H4.1|) in various quadrilateral do- 
mains have been obtained by the following method. A nonlinear algebraic coordinate 
transformation is first applied to map the quadrilateral domain on to the unit square. 
The equations are then discretized using centred second-order finite differences, resulting 
in a large linear system in block tridiagonal form. This is solved by a direct method and 
the total dissipation rate is evaluated, again to second-order accuracy, from a numerical 
integration based on either equation H2.9|l or equation p.lOfl (or their equivalents for 
the other differential equations). Typically equation (|2.10() gives a much more accurate 
numerical value for D, as it does not require a numerical differentiation of the solution. 

One of the simplest problems involving a wave attractor is to solve equation l|4.1|l in 
a tilted square domain (X = Y = 1, ti ~ t2 ~ t^ ~ t^ ^ 1) with uniform forcing 
(/ = 1). The dissipation rate for the case ti = 2 is plotted as a function of e in figure [3 
showing the anticipated convergence to a non-zero limit as e — > 0. Of course, a higher 
numerical resolution is required when e is reduced. FigureOlshows the spatial distribution 
of the dissipation rate for e — 0.01 and 0.003 as calculated with a numerical resolution 
of 1000^. As expected from the asymptotic solution, the dissipation is localized in the 
neighbourhood of the attractor and the width of the beam is proportional to e. 

To determine the expected limiting value of D from the asymptotic analytical solu- 
tion, some calculation is required. The focusing and forcing constants a and d are easily 
evaluated. A simple programme is then used to follow the propagation of rays in the 
quadrilateral domain and thereby to evaluate the ray-mapping and forcing functions 
Q^{9^) and A^{9^); these are shown in figure 0] 

Next the functional equations 

H^{e^{e^)) - H^{e+) ^ A^{e^) (5.1) 

must be solved, subject to the coupling condition on the Fourier coefficients 

H+^e^^'^-H-, n^O. (5.2) 
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Figure 2. Dissipation rate versus damping parameter for equation 14.11 1 in a tilted square 
domain {X = Y = 1, ti = 2) with uniform forcing (/ = 1). The dashed line indicates the 
limiting value expected on the basis of the asymptotic analysis. 



These colour figures are supplied as separate GIF files 



Figure 3. Spatial distribution of the dissipation rate for the problem referred to in figure |21 
The left and right panels are for e — 0.01 and 0.003, respectively. In each case a logarithmic 
colour scale over three orders of magnitude is used. 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Figure 4. The ray- mapping and forcing functions in the tilted square domain {X = Y = 1, 

ti — 2) with uniform forcing (/ = 1). 



Equations H5.1|l can be represented in the Fourier domain in the form 

oo 

^mriHn ^ H^n ^m: (5-3) 

n— — OO 

where 

M±„ = e-2'"-«%2nMe±(e±) jg±^ (5 4) 

Jo 

A,t=/ c-2™"'^*A±(0±)d0±. (5.5) 
Jo 

Equations (|5.3|l are semi-redundant because the symmetry properties Q^{Q^{9^)) = 0"^ 
and IS^(Q'^{9^)) = —/S.^{9^) ensure that the equations merely change sign when acted 
upon by the appropriate operator with matrix coefBcients It is most convenient for 
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Figure 5. Cross-section at height y — 0.3 of the j/-component of velocity, defined by 
Uy = —di/j/dx, for the problem referred to in figure|5|with e = 0.01. The upper panel shows the 
direct numerical solution, while the lower panel shows the sum of the asymptotic inner solutions 
for segments 2 and 4. The real and imaginary parts of Uy are given by the solid and dotted lines, 
respectively. 



numerical purposes to regard the unknown quantities as the 'large' Fourier coefficients 
(1 ^ n ^ A'^) and H~ {—N ^ n ^ —1) up to some finite truncation order N, 
substituting for the 'small' coefficients using equation l|5.2|l . The property = Smo 

means that the quantities cannot be determined and in any case are not required to 
calculate the asymptotic dissipation rate. In view of the semi-redundancy it is sufficient 
to solve only half of the truncated system of equations H5.3I) . e.g. the '-I-' equations for 
m > and the ' — ' equations for m < 0. The coefficients and are readily 

calculated to high accuracy from the computed functions 8^(6*^) and A^{9^). 

In practice it is found that, for a smooth forcing function /, the coefficients a„ typically 
decay rapidly with n so that oq makes by far the dominant contribution to D. This 
property is very convenient because ag is easily calculated and does not depend on the 
details of the global ray mapping. The dashed line in figure |2] indicates the expected 
asymptotic value oi D, to which terms other than oq contribute less than one per cent. 
Convergence to the same dissipation rate is found for the viscous problem defined by 
equation 1)8. l|l . 
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The asymptotic analytical solution is compared in detail with the direct numerical 
solution in Figure [3 which shows the y-component of the velocity along a cross-section 
through the tilted square domain at y — 0.3. In drawing the asymptotic solution along 
this section, the inner solutions corresponding to segments 2 and 4 are summed. No 
attempt is made to patch with the outer solution, which is much smaller and more 
difficult to determine properly. Even with this restriction the agreement is very good at 
e = 0.01. 

6. Discussion 

In this paper I have considered a prototypical forced linear wave equation featuring a 
wave attractor and weak viscous damping. This mathematical problem describes the long- 
term linear response of a typical bounded fluid to a periodic body force with a frequency 
within the range of inertial and gravity waves. Through asymptotic analysis confirmed 
by direct numerical calculations in an illustrative quadrilateral domain, I have shown 
that the forced disturbance is localized in the neighbourhood of the attractor and that 
the total dissipation rate is asymptotically independent of the viscosity. By considering 
a related equation with a non-viscous damping force, I have argued further that the 
asymptotic dissipation rate is independent of the nature, as well as the magnitude, of 
the weak damping mechanism. 

These findings have important consequences for physical problems such as that con- 
sidered in !j21and motivated by the experiments of Maas et al. (1997), in which a periodic 
body force excites internal gravity waves in a narrow tank of non-rectangular cross- 
section. In that problem the wave attractors exist in certain intervals of frequency that 
depend only on the geometry of the container. The long-term linear response of the fluid 
should give rise to a dissipation rate that is asymptotically independent of the viscosity 
and is a smoothly varying function of the forcing frequency within the interval associ- 
ated with a given attractor. This differs markedly from the case of an upright rectangular 
tank, in which the response of the fluid is governed by resonances with the regular invis- 
cid normal modes that can exist in that container, and the dissipation rate will be very 
small except where it is strongly amplified in the vicinity of the resonant frequencies of 
the lowest-order modes. As noted by Maas et al. (1997), wave attractors have a 'finite 
bandwidth' while the width of a normal-mode resonance is limited by viscosity. 

Ogilvie & Lin (2004) recently studied the linear response of a giant planet with an 
extended convective region to periodic tidal forcing. Inertial waves propagating in the 
annular convective region do not form regular eigenfunctions in the absence of viscosity. 
Numerical solutions indicated that, when the Ekman number is small, the disturbance 
tends to be localized on a web of rays. For intervals of frequency in which simple wave 
attractors exist, it was possible to identify their appearance in the solutions. In figure 10 of 
Ogilvie & Lin (2004) it can be seen that the dissipation rate appears to have converged 
to a value independent of the viscosity in the intervals where the dissipation rate is 
largest. It can be confirmed by tracing rays in the annulus that these intervals are those 
associated with the simplest and most powerfully focusing wave attractors. Although 
more complicated attractors exist outside these intervals, their focusing power is evidently 
not great enough to have led to a convergence of the dissipation rate at an Ekman number 
of 10~*. These numerical results are qualitatively consistent with the conclusions of the 
present paper, and in future work we will attempt to evaluate the asymptotic dissipation 
rate for tidal forcing in a spherical annulus directly using the methods described here. 

Forced inertial waves in a rotating incompressible fluid contained in a spherical annulus 
were computed previously by Tilgner (1999), who noted that wave attractors can be 
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detected in the solutions at small Ekman numbers but that other complicating features 
are also present. In his problem, which derives from earlier experimental studies, the 
forcing results from a sinusoidal modulation of the rotation rate of the boundaries, which 
is communicated to the fluid through viscous boundary layers. As the viscosity is reduced, 
the energy of the forced disturbance is also diminished. This behaviour is different from 
that found in the present paper with a distributed body force that is independent of 
viscosity, such as that resulting from tidal forcing. 

The apparent robustness of the dissipation rate is reassuring when it is considered that 
the true mechanism by which the inertial waves are damped is uncertain and difficult to 
model accurately. It is likely that the result is reasonably robust in other ways as well. 
For example, small corrugations of the boundary may not be so important for inertial 
waves as for acoustic waves, which undergo specular reflection. Inertial waves of a given 
frequency can propagate only in four different directions, and the angle of reflection will 
not be changed by a small corrugation. The wave attractors depend more on the gross 
shape of the container than on the precise contours of the boundary. 

It is also noteworthy that the 'inner' equations H3.30|l and (|4.4() studied in this paper 
are generic in the following sense. Provided that the inviscid wave equation is a hyperbolic 
second-order linear partial differential equation, its characteristics can be determined and 
any wave attractors identified. Curvilinear coordinates can then be introduced parallel 
and perpendicular to each segment of the attractor. If the equation contains a weak 
damping term resulting from a positive operator of either fourth or second order, then 
the coordinates can always be rescaled to obtain equation (|3.3Q|1 or equation (|4.4|l . re- 
spectively. This is true even if the inviscid equation has non-constant coefficients and 
the damping operator is anisotropic and depends on position, provided that there are no 
singularities such as corotation resonances or critical latitudes. 

Further work is required to understand the transition from normal modes to wave 
attractors as the geometry of the container is varied. Ogilvie & Lin (2004) found that 
some 'memory' of the normal modes of inertial waves in a full sphere was retained when 
a small solid core was introduced. Even though no inviscid normal modes are believed 
to exist in a spherical annulus, the dissipation rate is still enhanced in the vicinity of 
the eigenfrequencies of the lowest-order modes of the full sphere. An imperfect resonance 
is possible with such modes, which have a length-scale larger than the radius of the 
core and are not greatly affected by it. At the same time, introduction of the solid core 
allows a much richer response at other frequencies, associated with wave attractors. It is 
possible that in such problems the attractors achieve complete dominance only at Ekman 
numbers yet smaller than those achieved in recent numerical studies. 

A remarkable theory of axisymmetric inertial waves in an incompressible rotating fluid 
contained in a spherical annulus was developed by Rieutord et al. (2002), who considered 
normal modes that are regularized and damped by viscosity. Their problem is similar to 
that considered in the present paper, but there are also important differences. Neglecting 
some of the effects of curvature so as to make the equations somewhat simpler, they 
found asymptotic solutions for normal modes in the limit of small Ekman number, in 
excellent agreement with their numerical calculations. The modes are localized near a 
wave attractor and satisfy an eigenvalue problem analogous to that for the wavefunction 
of a quantum-mechanical harmonic oscillator. The width of the shear layers in these 
modes scales with v^^^, although the authors noted that ly^^^ layers also play a role when 
the problem is studied in full spherical geometry. In contrast, the solutions described in 
l^lhave only a v^^^ shear layer and satisfy a different inner equation that has no strictly 
localized solutions but is forced from the outside. The reason for the different scalings 
is that the normal modes considered by Rieutord et al. (2002) have frequencies that 
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are asymptotically close to the edge of the band within which the attractor exists, so 
that the attractor is very weak and its focusing power a is asymptotically close to 1. In 
contrast, the forced solutions in the present paper are constructed for the more interesting 
frequencies that lie properly within the bandwidth of the attractor, so that a is not very 
close to 1. It is likely that a different asymptotic theory of forced oscillations could be 
constructed for values of a very close to 1, which would more closely resemble the analysis 
of Rieutord et al. (2002). However, it still remains to be demonstrated in detail how the 
analysis in the present paper applies to forced oscillations in spherical geometry. 

I thank Doug Lin for introducing me to this problem and for ongoing discussions, as 
well as for hospitality at UC Santa Cruz where some of this work was carried out. I also 
thank Raincr HoUcrbach, John Papaloizou, Michel Rieutord and Yanqin Wu for helpful 
discussions, and the anonymous referees for their useful suggestions. I acknowledge the 
support of the Royal Society through a University Research Fellowship. 



Appendix A. Character of the inviscid linearized equations 

The equations governing the dynamics of an ideal fluid with negligible self-gravity may 
be written 

^ = -lvp-V$, (Al) 
gf = -pV.«, (A2) 

^ = -7pV.'u, (A3) 

where $ is the external gravitational potential and the notation is standard. Let (r, </>, z) 
be cylindrical polar coordinates and consider a basic state consisting of an axisymmetric 
fluid body with angular velocity Q{r, z) and no meridional flow. For small perturbations 
of the form 

Re [u'(r,^)e-'"*+''"'^] , (A 4) 

etc., the linearized equations read 

I f , d , d \ , imp' 
iL— + u'— r^Cl) = — 



-™;=-i|::+4l^, (A 7) 



p dz (p- dz^ 



-— r< + ^ + 

r or r oz _ 

Id. imu'. du' 
-— r< + - + ^ 

r or r oz 



(A 8) 



(A 9) 



dr ^ dz 

where Cj = u — mfl is the Doppler-shifted wave frequency. Eliminate it^ and p' to obtain 

ic,^ - AK - B< = f ^ - f^^) + 2n^, (A 10) 

p \or or jp J pr 
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-CK + iu^^-DW^^-'-^f'^-^^], (All) 
217 9 ^ 2oN I dp f I dp 1 dp 



where 



r or p or \jp or p Or J 

r az p or \^p oz poz J 



p dz \ jp dr p dr 

I dp f 1 dp 1 dp 
p dz \"fp dz p dz 
Then ehminate u'^ and u'^ to obtain 

.d)^ + {B + C)^ + (c2;2 _ + (terms in p' and Vp') = 0. (A 16) 

This equation is hyperbohc when 

'i{Lo'^ - A){Cj'^ - D) < {B + Cf , (A 17) 

i.e. for frequencies such thatf 

{A + Df - [{B + Cf + {A- Df] < 2^2 < (A + i?)^ + [{B + Cf + {A- Of] 



(A 18) 



Appendix B. Vortical effective forcing 

Ogilvie & Lin (2004) studied the hnear response of a rotating giant planet (or star) to 
tidal forcing in circumstances in which the forcing frequency is comparable to the angular 
velocity Q. of the body but small compared to its dynamical frequency [GM / B?)^^"^ . They 
used a consistent asymptotic expansion to simplify the equations satisfied by the tidal 
disturbance in both convective and radiative regions. The tidal force, which derives from a 
perturbing gravitational potential, is taken up by a quasi-hydrostatic adjustment (a tidal 
bulge) known as the equilibrium tide. Because the tidal frequency is non-zero, the motion 
of the bulge gives rise to a non-zero velocity field associated with the equilibrium tide, 
which can be derived analytically. The existence of inertial forces associated with this 
flow mean that the quasi-hydrostatic bulge does not quite satisfy the equation of motion 
and a residual force appears. 

In quantitative terms, the total velocity perturbation u in an adiabatically stratified 
(i.e. convective) region satisfies the equations 

- iCju + 2Q, X u = -VW, (Bl) 

M + v.(pM)=0, (B2) 

where Cj = lu — mil is the Doppler-shifted frequency as in Appendix^ W is the pressure 
perturbation divided by the density, plus the total gravitational potential perturbation, 
and is the density perturbation associated with the equilibrium tide (any additional 

f The H0iland criteria for stability with respect to axisymmetric perturbations are A + D > 
and AD > BC (e.g. Tassoul 1978). For a H0iIand-stabIe basic state the hyperbolic range of 
frequencies may or may not include a) = 0. 



Wave attractors and tidal disturbances 25 

density perturbation being much smaller). Here it is assumed for simplicity that the fluid 
is uniformly rotating and the viscosity is neglected. Now since 

^ + V.(pue)=0, (B3) 

the residual 'dynamical tide' Ud = u — Ue is found to satisfy the equations 

- iu)Ud + 2n X Ud = -VW + a, (B4) 

V-ipud)=0, (B5) 

where 

a = iu)Ue — 20 X Ug. (B 6) 

The dynamical tide therefore obeys the anclastic approximation, in which the modified 
pressure perturbation W must adjust to satisfy the constraint V • (ptta) = 0, and is 
driven by an effective body force a per unit mass. In general this effective force is vortical, 
because V x a 7^ 0. 
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